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MATHEMATIC MODELING OF THE EARTH'S SURFACE AND THE 
PROCESS OF REMOTE SENSING 


B. M. Balter 


It is shown that real data from remote sensing of 
the earth from outer space are not best suited to the 
search for optimal procedures with which to process such 
data. To work out the procedures we propose to use 
data synthesized with the help of mathematical modeling. 
Subject to the condition that these data are similar to 
the real data, the processing procedure worked out for 
them will also work well in practice. A criterion for 
similarity to reality is formulated. It is shown that 
a simple method for staisfying it exists. The basic 
principles for constructing methods for modeling the 
data from remote sensing are recommedded. A concrete 
method is worked out for modeling a complete cycle of 
radiation transformations in remote sensing. The crux 
of this method is a suitable procedure for modeling a 
multi-dimensional random field. We describe a computer 
program which realizes the proposed method. Some re- 
sults from calculations are presented which show that 
the method satisfies the requirements imposed on it. 

The required statistical characteristics of the terres- 
trial surface and the radiation transformations are 
assigned to the program anu it produces a frame for 
remote sensing which is ready for processing. 


INTRODUCTION 


The data from remote sensing are not the best raw material for 
experiments . In remote sensing (RS) the sole observable quantity, 
from which information can be extracted concerning the condition 
of the terrestrial suface, is the reflectance of the surface. 

(Here we do not consider natural radiation and polarization, since 
taking them into account changes nothing in principle.) As is 
well-known, reflectance depends on a multiplicity of factors: 
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soil, the vegetation phase of flora, humidity, etc. These factors 
are extremely variable, are often difficult to measure directly, 
and are completely uncontrollable. It proves to be impossible to 
obtain two images of parts of the surface for which these factors 
are identical, or in which only one of them is different. Cf. 

[1, 2, 3]. Therefore we cannot study the influence of individual 
factors on the quality with which the frame for some precedure is 
processed, and the investigator turns out to be in the position 
of a person who attempts to evaluate the voice of a singer by 
listening to his singing in a chorus. Due to the difference in 
the original data it is impossible to compare the results of pro- 
cessing obtained by differing groups of investigators. Cf. [1, 4], 
Because of the large number of factors it has not been possible to 
find the reasons for the difference in the behavior of the processing 
procedures for various regions [5, 6, 7]. 

Note that in the exact sciences, for example, in physics, the 
ability to repeat, and the strictness with which conditions are 
controlled, are unalterable requirements for an experiment. If 
a physicist were asked to work in a situation in which the charac- /4 
teristics of the patterns (frequency) and the conditions of an 
experiment (temperature, magnetic field) do not depend on him, 
cannot be reproduced at will, and are rapidly changing, he would 
doubtless regard such conditions as unacceptable. This is precisely 
the situation when we have occasion to work with remote sensing. 

Certain methods of processing: [5, 8, 9, 10, 11] are used for 
training in the procedures for recognizing a priori information 
from regions which are selected as standard for an entire frame. 

This information is adequate only for a certain neighborhood of 
the standard regions. In such methods it is absolutely necessary 
to estimate the influence'- of the degree of adequacy and the complete- 
ness of the a priori information on the quality of the work [1]. 

But in order to do this, it is necessary to have a priori informa- 
tion from the entire frame, which is impossible in practice. 
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The traditional methods for working out the procedures in RS 
are reminescent of an attempt to learn English or French not from 
textbooks, but from the works of Joyce or Proust. Even if it is 
necessary to read precisely these authors, it is not necessary to 
learn the language from these works. And in RS, although it is 
necessary to process the real data, it is better to learn processing 
from some type of textbook. In the present paper we also propose 
a method for creating such textbooks. We refer to mathematical 
modeling of data, so that the results resemble the real data in 
their basic features, but have controlling mathematical characteris- 
tics. After checking out the procedures on synthesized data it 
will be possible to carry out a "precise adjustment" to the real 
data. 


Modeling of data. Advantages. A Quality criterion . Replacing 
natural data by artificial data makes it possible to change some 
characteristics of the frame being processed without touching 
the others. Moreover, precise information about the entire frame 
is available to the experimenter. However, we must ensure a 
connection between the results of testing some processing procedure 
and reality. Indeed, in typical models of signals and their re- 
presentations, the data is obtained in analytical form (e.g., 
additive Gaussian noise is chosen) and analytical evaluations of 
quality are obtained, but they do not guarantee quality under 
real conditions. See [1, 12, .4], It is clear that a criterion is 
necessary for the quality of the modeling which ensures conver- 
gence of the process of optimization for the procedure to a proce- 
dure that is optimal for the real data. Such a criterion is for- 
mulated in the following section. 

The author knows of only one paper on modeling in RS. This 
is [13] in which the change in the radar image of a locality is 
simulated when the observation angle and the form of the projec- 
tion are changed. The basic purpose of [13] was to detect parts 
of the locality which were concealed from the beam from a radar 
station. The reflection coefficient at each point is assumed given. 
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Therefore [13] corresponds to one unit of the given method: 
modeling the process of scanning and conversion of radiation into 
a signal. The possibilities for applying modeling . The various 
applications of artificial data can be divided into A groups. In 
the first place, it is possible to perserve unchanged the modeled 
frame and to select the optimal paramters of the photography 
(height of the flight, size of the frame, settling on the locality, 
the height of the sun, etc.). This might also be done with real 
data but for the Impossibility of keeping the landscape unchanged 
while it is being photographed several times in different ways. 

Secondly, It is possible to change the properties of the noise 
(static, the noise of the receiver) and to determine the domain 
of effectiveness of the processing procedure. In the third place 
and this is the most interesting group of applications, it is 
possible to change the set of objects on the terrestrial surface 
and their characteristics: humidity, soil type, contamination, etc., 
and for each set of characteristics to compare the various process- 
ing procedures. Then we shall clarify to which changes in the 
framework under observation the quality of the processing is most 
sensitive. Moreover, we shall be able to decompose the region 
of changed characteristics into subregions, inside of each of 
which some precedure is optimal. Knowledge of such a decomposition 
makes it possible in practice to choose the most suitable process- 
ing procedure. In the fourth place, since the entire a priori /6 

information is at our disposal, we can change its amount from 
0 % to 100/S and can clarify which minimum is necessary to ensure an 
acceptable quality of work. This makes it possible to determine 
the structure of the ground measurements and the data bank. 

In the general case it is possible to investigate some proce- 
dure with respect to all four directions. This makes it possible 
to check it out in the laboratory, so that in the case of real work 
it is only necessary to do a little polishing, which can be carried 
out by the method of self-instruction. Clearly, it is impossible 
to undertake the investigations enumerated with respect to real 
data. Up to now such attempts have only lead to the accumulation 
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of piles of useless information. However, as is clear from what 
has been said, the use of the methods of modeling might place the 
processing of RS data on a more or less rigorous basis. 

The method proposed in the present paper makes it possible 
to carry out all the indicated investigations. 

1. THE PRINCIPLES OF MODELING IN RS 

The transition from recognition to measurement . For the 
time being,- the processing of data in RS is limited basically to 
recognition of ground objects and the construction of thematic 
maps. The search for effective recognition procedures is incomplete 
(to a significant degree due to the problems enumerated above), but 
the role of RS is not restricted to recognition. The basic purpose 
of RS Is the measurement of the most important characteristics of 
objects (productivity, pollution, etc.), since this is most impor- 
tant economically. The measurement must be sufficiently precise 
so as to be competitive with the traditional methods for obtaining 
the same information. There are still no methods for evaluating 
the quality of measurement procedures, not the procedures themselves, 
but assuming that there is a future for them, we have oriented 
the modeling twward them. The ideology of the parameters which 
is proposed in the following section serves as a basis for modeling 
the characteristics being measured. It seems that in this or the /7 

other form is should inevitably appear, even in the measurement 
procedures themselves. The description of objects with the help of 
the parameters . Presently the frames in RS are simply described 
by the luminosity at each point, i.e., directly by an observable 
quantity. Such a phenomenological description is characteristic 
of all sciences in their early stages. But it would be practically 
impossible to study, e.g., hydrodynamics or the molecular theory 
of heat if we were only able to describe motion phenomenologically, 
going into all the details of the complicated changes In positions. 
Fortunately the characteristics of bodies which determine their 
motion uniquely (mass, velocity, acceleration, external force) are 
known . 
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A quite similar conversion is also possible for RS. It is 
well-known that humidity, the phase of vegetation and the soil 
strongly influence the reflectivity of agricultural crops, and 
that depth, pollution, etc., affects the reflectivity of reservoirs, 

[14, 15]. We may assume that the specification of a sufficiently 
large set of such characteristics ("parameters") with great pre- 
cision determines its reflectivity. (For brevity we shall call 
reflectivity albedo , although it depends on wave length and the 
observation angle.) Once having measured the dependence of the 
albedo of objects on their parameters, it is then possible by 
describing the frame in the language of the parameters to recover 
its observable appearance. Of course, for different objects, the 
parameters are different. 

It is necessary to impose two requirements on the parameters. 

First, they must uniquely determine (within the limits of the speci- 
fied precision) the albedo of an object. Otherwise, the description 
of the frame with their help will not be complete. Second, they 
must be intrinsic, i.e., they must correspond to the characteris- 
tics generally accepted in geography, the agricultural sciences, 
etc. Then it will be possible to use the information about the 
parameters accumulated by these sciences. 

The quantities which RS serves to measure: productivity, pol- 
lution, etc., affect the albedo (otherwise it would be impossible /8 

to measure them), and this means that they fall into the category 
of parameters. A parametric representation makes it possible to 
introduce these quantities into the processing in explicit form 
and to describe the framework as a system of stochastic equations 
relative to the parameters. The quality of a processing procedure 
is naturally characterized by the precision with which the para- 
meters are measured. Note that the paameters are spectrally 
invariant, and therefore the values of the intensity in different 
channels are easily associated in terms of them. 

It seems that a parametric description arises naturally in 
RS. The proposed method is entirely based on such an ideology. 
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The statistical approach . In RS the statistical treatment of 
data now predominates. Accordingly we shall assume that the treat- 
ment procedures produce the statistical characteristics of the para- 
meters being measured, e.g., their histograms of the distribution 
of probabilities. Experiments with the method of modeling consits 
in changing the desired statistical characteristics of the para- 
meters and transformations of radiation (the "input data" of the 
method) and recording the corresponding changes in the output of 
the processing procedures. The criteria for the similarity of 
the modeled data to reality also have to be formulated in statisti- 
cal terms. 

Specific criteria for the quality of the modeling . We shall 
call the precision with which the statistical characteristics of 
the parameters being measured are recovered with the aid of the 
processing procedure the quality of the procedure. To model a frame 
the necessrry input information is: the desired statistical charac- 
teristics of the parameters, etc. The synthesized data correspond- 
ing to some real region of the surface are data obtained by modeling 
in accordance with the input information gathered in this region. 
Naturally the modeling is considered successful if the quality of 
the processing procedure is considered identical for both the real 
and the corresponding modeled data. Now we shall formulate the 
criterion for the quality of the modeling . 

Let U be a region of the surface. Let us photograph, it from /9 
outer space and process the photograph by some procedure T. The 
procedure will provide histograms of the probability distributions 
of the parameters i: = T(U). Let us collect on U the 'input data 

for the modeling U and with respect to them and with the help of 

O 

the method M, let us model the frame W = M(U„). Processing this 

° i 

frame with the same procedure T will give some histograms P 5 . Now 

i i ^ 

let us compare P^ and P 2 and ascertain whether they reflect the 

same (i.e., the "true") distribution of the probabilities. For 

this purpose numerous criteria have been developed (the sign tests 

of Wilcoxson, Kolmogorov, etc.). When these criteria are satisfied 
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for all 1 or for some function of particular criteria, we shall 
consider that the quality criterion for the modeling is satisfied 
for the procedure T and the frame of U. 


Two defects remain: it is necessary to verify the criterion 
anew for each new procedure and for all input data, while drawing 
on the corresponding real frames (and where do they come from?). 

It is desirable to verify the method only once, even if this in- 
volves a complex method, so that this verification ensures satis- 
faction of the quality criterion for sufficiently broad sets of 
processing procedures and input data. We shall show that this is 
possible . 

We shall assume that the description of the region of the ter- 
restrial surface U and the transformations of the radiation include 
a frame observable from outer space and a priori informatinn about 
the parameters of ground objects and the atmosphere, and also the 
laws for transforming the parameters into the observed radiation 
intensity. This description (which possesses some redundancy) 
generates a space of frames A. It is possible, for example, to 
specify the parameters and the intensity values of the radiation in 
all the spectral channels, for each point of the frame A is compact, 
i.e., the number of points is finite and the values of the radia- 
tion intensity and the parameters are bounded. 

Let us introduce a metric d into A. As is well-known, any 
compact set can be covered by an e-network with a finite number of 
points U i such that for any (Jef A there exists a U.^ so that 

c[(U ii V)&£ (i- 

Note that the criteria for the agreement of the distributions, 
which are based on quadratic functions, possess metric properties. 

Since continuous transformations do not affect compactness, 
we iudy choose any description of the frame from a broad class which 
preserves compactness. We shall select the description which has 
the most convenient e-network. We shall denote the set { } by 


/10 


8 



S and we shall use It as Input data for modeling. We shall describe 
the method of modeling by the operator M which maps S onto A. The 
operator M is stochastic, but here we shall require that each of its 
realizations satisfy the quality criterion. The result of process- 
ing is the subset S, and to the procedure of data processing corres- 
ponds the operator T, from A in. S. An illustration of the struc- 
ture described is given in Figure 1 (p. 27). 

The criterion introduced above for the frames U and W with 
the significance level a is described by the condition 

1 [f (U), T(W)] £ «< 'r- : V'f= W(Vs);. : d - 2 


U is the projection of U onto S, i.e., input data collected 
s 

on the frame U and analogously for W . Let d(W , U ) = 3. Let 

s s s 

us impose on the procedure T the Lipschitz condition relative to 
A with constant K (which can be evaluated analytically), uniform 
on the class of the procedures under consideration: 

' d(r(J,) l ' r(u 2 ))< Kl(u,,(A). u.b) 

This means that the procedure does not yield strongly differing 
results for more or less close frames. Then 


l[T(U);?(W)]< Kl(uy)< K[i(ii,u s )4K^)fl(W,y)]^K (Z£+f), 


(1.4) 


It is clear that if a is given, then provided that e is chosen 
sufficiently small, i.e., the volume of input data is sufficiently 
large, so that on the class of the procedures used 2Ke < a, we 
obtain a sufficient condition for the satisfaction of the quality 
criterion in the form K (2.S+ ft) ^ . That is, for a given class 

of procedures and a threshold a we select a volume of input data, 
after which the quality criterion reduces to the condition 

, (1.5) 

in other words, to testing the precision of the reproduction of 
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the input data in the modeled framework. This test is purely mathe 
matical and does not require real data. The latter are used only 
in the selection of s. 


Thus we have shown that it is possible to simplify the criter- 
ion. This is not, a rigorous proof, but rather a proof scheme. 

Some practical questions . As a priori input data for modeling 
we have selected the following. 


1. The marginal density of the probability of the parameters 

for ground objects P^ (X ) where a^ is the ith parameter of the 
jth object. The coefficients of the correlation between the dif- 
ferent parameters at one point of the surface i and the homo- 
geneous, isc tropic functions of the correlation between the values 
of a single parameter at different points of the surface (x^ y^ 
and(x 2 ,y 2 )= ^ ' 

2. The analogous densities of the probability and the correla- 
tions f or the parameters of the atmosphere: (X) t Kz and 

pj' ()h~ M * Here K indexes the parameters of 

the atmospnere in the spectral channel l . 


3. It is necessary to measure, and once for all build into 

the method, the functions which describe the dependence of the spec- 
tral reflection factor and the indicatrices of scattering on the 
parameters for all ground objects, and likewise for the dependence 
of scattering and absorption in the atmosphere on its parameters. 

4. The remaining input data (the characteristics of scanning, 
of the receivers of radiation, the general location of objects) 
are less important. 

These input data were chosen because they are already being 
measured in different sciences, and to measure any others is diffi- 
cult in practice. Starting from an experiment it was assumed that 
such data form a sufficiently dense e-network on the set of 
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frameworks. Of course, the quantity e must be determined experimen- 
tally, and it is limited by the attainable quality of the modeling. 

Note that the criterion for 0 is only a sufficient condition. 
Within its scope, it is possible to improve the agreement with 
reality, proceeding from empiricial considerations. We shall assume 
that the most important objective is to model precisely the densities 

r ; ^ 

: t -7 and 

a i ' - ■ . — ■ — ^ _ _ * — j / *- 

f I y j » and then Pj^ (x) for the atmosphere, 

and finally, the correlations of the parameters of the atmosphere. 

The data in group 3 are assumed to be precise (this depends only 
on the quality of the ground measurements). This ranking in terms 
of importance proves to be essential. 

The proposed method possesses a definite generality, but 
basically it is oriented toward photography of agricultural crops. 
This is the most economically effective application of RS. In geo- 
graphical problems, for example, there exist geometrical structures 
to which statistical ideology is not applicable. 

In the following two sections we describe a concrete realization 
of the method of modeling based on the concepts discussed above. 

2. A DESCRIPTION OF THE METHOD 


The basic steps . The method consists of the following basic 
steps . 


1. Modeling a plausible arrangement of the objects or reading 
it out from an external lay-out. 

2. Modeling the values of the parameters at each point. This 
is a random field with characteristics which are given in the input 
data. 


3. Transformation of the parameters into the albedo in accor- 
dance with given functions. Modeling the field of the radiation 
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which is incident on the earth, and transforming it into the field 
of the escaping radiation with regard to the albedo. 


4. Modeling the field of the atmospheric parameters. The 
corresponding transformation at each point of the frame. 

5. The imitation of scanning and the transformation of radia- 
tion into the RS singal. 

A more detailed block diagram of the method is presented in 
Figure 2 (p. 28). Let us discuss in the order of importance (steps 
2, 3, 5, 1) the basic problems and the solutions. 


Modeling the parameters: requirements . The random field of the 
parameters a^ 1 where j = 1, ..., M indexes the objects, and i = 
1,..., N . , the parameters is multidimensional. We shall require 
that the field have given marginal distributions of the probabili- 
ties PoA 00 and have given correlations /*• A and 


• ■ — ' — n (<2 

j'V i } r “ X 3 * Here restrictions must be imposed 


on 


the form of the distributions. 


This requirement must also be satisfied in case it is neces- 
sary to confer previously assigned values on certain parameters 
at certain points. Indeed, such a parameter as the height of a 
locality can be completely fixed in advance. Moreover, we may fix 
a part of the parameters, and alter the other part. Finally, when 
investigating the influence of standard information it is possible 
to fix the parameters on the standard region and alter them on the 
remaining territory. 


In practice the field of the parameters is not stationary, since 
their statistical characteristics have drift along the framework. 

This is important since drift limits the maximum size of a region 
which can be processed in one reception and determines the density 
of the network of points in which it is necessary to collect a 
priori ground information for the processing procedures. Therefore 
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it is necessary to model a nonstationary field also. 


Finally, the quickness in the response of this step is also im- 
portant. It is potentially the most voluminous, since at each point 
it is necessary to compute the correlations of the parameters with 
respect to its neighborhood (possibly, quite large). 

In toto we make noticeably more stringent demands on a random 
field than usual [16, 17]. We propose that the modeling procedure 
satisfy them. 


Modeling the parameters: realization. If N, is the number of 

1 J 

parameters a J 1 for an object j, we shall select N = max(Nj) and 
we shall model an N-dimensional random process. For points with 
different objects its components correspond in general to different 
parameters, but this complication is compensated for by the number 

of components from £ /» . to max(N. ) which provides great computa- 

J J 

tional economy. In order that it may be possible to work effec- 
tively with arbitrary J. (X') let us perform the following trans- 
formation. In place of the desired parameters let us model the 
field of the normally (0, 1) distributed random quantities 
Let us denote their probability distribution functions by Q wJ . f X , 

and the correlations by j)^ ^ and ^ • 


Then we note that the transformation 


(Ju J £ = f 1 ^ '[6^ (f £)] where X . (2.1) 


yields precisely the desired distribution of the probabilities 
P r j. (X; . The correlations P are found from 

U L 

# % 

r- [aUl M [>£,1 

~ V^L'a^]Afa{J 

where M and 2 are the mathematical expectation and dispersion 
operators, respectively, and 
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(2.3) 


co , , I 


Thus we obtain a convenient modeling scheme. 

1. Given to find and $ [jCi. ^ ] . 

2. Given f* to solve the system 2.2, 2.3 for p. 

3. To model the field with Qaussian (0, 1) Joint portability 
density and given correlations p. 


4. To carry out at each point of the field the inverse of 
the transformation 2.1. 


The laborious steps 2.2 and 2.3 are performed only once , 
independent of the number of points in the frame, but the modeling 
at each point of the Gaussian quantities is carried out quite 
rapidly. For the p (a) there are standard approximations and 


typical forms, for whicii st^p 4 (i.e., 2.1) also takes up an accept- 
able amount of time. 


Generally speaking, equations 2.2 and 2.3 may not have solu- 
tions. But if the distributions are not too far removed 

from normal distributions, or if the correlations are not too 
close to unity, this risk does not threaten. Cf. [17]. 

In such a scheme it is easy to take account of a nonstationary 
condition. For this purpose transformation 2.1 is carried out 

• .-t 1 ^ 

with a new J • , taking account of drift. It is true that p 

* u* t 

was computed only for a single f , so that the obtained 

by the inverse transformation will have correlations somewhat dif- 
ferent from the given one. But, as is clear from 2.2 and 2.3, the 
error has the same order of smallness as the nonstandardness of 
; thus it can be reduced to an arbitrary value of the frame 
processing "in the squares". Moreover, we have arranged the 


characteristics with respect to importance and we shall give pre- 
ference to precise reproduction of the drift of the distributions 
over precise reproduction of the correlations. 


There remains the problem of fixing certain values of a^ ^ . We 
shall solve it by a method which in general significantly increases 
the flexibility of the modeling for a random field. We propose to 
replace the traditional raffle of all N components of the process 
at once at each point by a successive modeling of conditional 
distributions. In fact, 




where S denotes the values in the neighborhood of the 

given point. If certain values ^ £ are fixed, we must simply by- 

pass them when modeling, and the given form f ^ a/) will 

be maintained automatically. In the neighborhood of n/ c only 
those points are included which contain the very same object as 
the neighborhood under consideration. Note that some hierarchy is 
inherent in the parameters (the topography depends on nothing, the 
soil depends only on the topography, etc.) which the successive 
modeling approaches. 


The key feature in the sequential scheme is the possibility 
of calculating in advance the form of the dependence of V] 

. j i j . j * \ I K j // 

V * . , S ) on S and ^ A . This is impossible for arbitrary 

distributions. But in the present case the joint distribution, 

and therefore all the conditional distributions t J are normal. 

5 t 


Then, in accordance with [18], for the jointly normal (0, 1) 






is normal with 


where 


A 


M M - - £ 

is the determinant. 


Arc 


A- 


/\k>. 


(2.5) 


and /\ k t , the cofactor of the element 
/v. k £ in the matrix of the second moments of the joint distribution 
That is, in the distributions 2.4 the dispersion is known. 
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and the calculation of the mathematical expectation reduces to a 

• • 

summation with given weights with respect to . and 

J fj •• i 5 

S, which is extremely simple. Therefore the combination of the 
transformation 2.1 with a sequential scheme is successful. 


Next between stages 2 and 3 in the modeling of a random 
field it is necessary to calculate /\^i . For £ J the columns of 

J • • 

II All correspond a) to the quantity ^ ^ itself; b) to 

at the very same point; c) to ^ ^ in the neighborhood of S; d) 

v, ; A 

to $ j .•••*$/- d in S. The computation of ''Kl reduces to the inver- 

sion of U f\ |J ; A,U” li^ll • The inversion of a large 

matrix is quite laborious. But it is logical to assume that the 
parameters of the group b) are connected to c) only indirectly in 
terms of a). In the general case, for three groups of quantities 
Y i, Y 2 , Y 3 the condition for Y-j^ to depend on Y^ in terms of Y 2 
only is written as 


P ( Y<, Y, | Yi) - P(Yi I V-V f ( Y* I Yi) 


( 2 . 6 ) 


Then according to the theorem of Bayes 


r , v iy y\ p(^>ir»L - umivw' jwiW’Ww fliPhA P 

pc v-v,) ■ -wy 


7 ) 


Thus d) can be discarded, and .>Ki can be calculated with respect 
to the submatrix consisting of a), b) and c). Further simplifications 
are not possible. 


Thus a method has been constructed which meets the demands of 
the preceeding section. Intrinsic modeling reduces to the computa- 
tion of the linear form (2.5). 


Modeling the atmosphere: problems ; 1 . As is well-known, there /17 
is as yet no unique model for the transformation of radiation in the 
atmosphere. This means that it is necessary to construct a model 
such that a change in the model causes only a change in the coeffi- 
cients which ocour in the method. 
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2. Due to the limitations on computing time, it is impossible 
to take into account all the components of the atmosphere which 
affect radiation. 


3. The behavior of radiation depends on the distribution of the 
atmospheric components along a ray. In practice this information 

is rarely available. 

4. From each point the scattered radiation is spread over a 
large neighborhood in view of the great length of the path of the 
ray and the multiplicity of the scattering. Computer time does 
not permit calculation of this entire neighborhood. 


Modeling the atmosphere: realization . In the statistical 
approach these problems are significantly moderated. The task 
reduces to computing the change in the distribution of the proba- 
bility and in correlations of the radiation during passage through 
the atmosphere when the distributions of the probability and the 
correlations for the parameters of the atmosphere are given. The 
following computational method is proposed which is capable of 
accommodating various concrete models of the atmosphere. 


The parameters of the atmosphere are the density of the air 
the 

diffusers 

■ *, 

m 

23 


d and the concentration of the "effective" absorbers q and the 

0 

p^ of radiation in each spectral channel £ i * 

tjsiixvi a n'ujmJ* 


ji 

>_ 

c-1 JU 




^ 3 S 33 a 


( 2 . 8 ) 


Here t 1 and u^^ are the concentrations of each of the n diffusers 
and the m absorbent components of the atmosphere; ls the 

Ith spectral channel; A the wave length; g and f the coefficients 
of the dispersion and absorption. The atmosphere is assumed to be 
plane and infinitely thin, so that Pc, and are the mass of a . 
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unit area. Generally speaking, this approximation is extremely 

crude, and the statistical characteristics are slightly distorted, 

if, as necessary, we select the functions expressing the dependence /l8 

of the scattering and absorption on X .1/ and d c . These 

functions and 6 c ! (a,4'£) must be determined from a concrete 

model. In the first approximation we may set 


« Jr At[£ 

ir- iTi m[| 4 


( 2 . 9 ) 


The integration is performed along a column of atmosphere of height 
L. The random variables are t, u and d. For a thin layer LZ the 
intensity of the previous radiation I^LZ) in the channel ( is 


1 


( c ) Xiioii- P c h^f : J 12$ 


Analogously, if is the radiation at the "entrance" to the 

atmosphere, then in our model we shall define the radiation at the 
"entrance" as 


1^ (*, a)- I t Mil M ^ (fe'A -i(vO l] k (k'A i, A' j )> 


(2.10) 


/ ^ ^ ^ p (x y) and ^ are modeled as the random fields • 


The scattered radiation is computed as follows. The 

distribution of the probability , the correlations 

ii t _ j * C' 

' L\ and i*£ j\ V J are com P uted with respect to 

,j / X .*) and A • Then the field of the random variable J L (*,b) , 

with characteristics which depend on (X, and H [*X* (x ( b)]] * i s 
raffled: 


Pj, (k,o l\ c<.,) W] ' M ft (*■ 

i\c<* = ft c ,,)] 


( 2 . 11 ) 


Here the absence of local characteristics is connected with the 
assumption concerning the great length of the path and the multiple 
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scattering of the radiation. The total radiation incident on a 
point of the earth (x, y) or incident from the earth onto a satellite 
is 


I." (x,s)= T.l M - 3 t A [ </ M.-m Piftsrtl} 


( 2 . 12 ) 


Here is the indicatrix of atmospheric scattering which 

depends on the angle 4> between the line of sight and the direction 

of the average radiation current. /19 


Thus modelling the atmosphere reduces to the following steps. 


1. Raffling the field of the characteristics of the atmospheric 
parameters with respect to these characteristics. 


2. Modeling the field of the scattered radiation by means 
of formulas (2.11). 


3. Calculating the input radiation I£(x, y) by means of for- 
mulas (2.10), (2.12) and the input radiation (I^(x, y) 

By suitably choosing L t - [L.t'i) and , it is possible 

to obtain any needed distribution of the probabilities of the output 
radiation. As shown above, we attach less importance to the corre- 
lations. The given model is plausible, but it needs to be tested 
experimentally . 


Transformation of the parameters into albedo and the escaping 
radiation . Here we take into account the indicatrix of the scatter- 
ing, the illumination of the soil by radiation passing through 
vegetation, the presence of shadows, etc. The transformation is 
described by the following formulas. 


iiM » (2.13) 

where > ;k u =Je 5 

h‘(u LAjx.i) 5 • 

- pdi fini is the P r °J ective covering of the soil by the object j; 

IV ~ "iosr 
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fid} is the density 0 / the planting which has a given probability 
distribution; is the projective covering when **I ; 

^ is the inclination angle of the locality; are the Para- 

meters of the object j; £, n are angles in the - plane connected 
with the sun (cf. [14]); £ 0 and n 0 are £ and n for the satellite to 
point (x, y) line of sight; (j r^j is the angle between ^ 

and the direction (x, y) to tne^sun; is the indicatrix of scat- 
tering for the object, and )C that for the atmosphere, 4 / is 
their convolution, which plays the role of the indicatrix of the 
object with respect to the scattered light; is the albedo 

of the object; the albedo of the soil. The functions F£, 

G£, are determined experimentally and are once and for all incor- 
porated into the method. I^(x, y) is the direct, and (X t y) 

the scattered solar radiation incident on the earth; I^(x, y) 
is the total radiation escaping from a point in the direction of 
the satellite; a = 1 in the sun, and a = 0 in shadow. Thus the 
transformation consists of the following steps. 

1. Calculation of the field of the incident light and the 
indicatrix of the atmosphere. 


2. Calculation of the object’s albedo and the albedo of the 

soil by means of the given functions F£ and and the values of 

parameters n J , n J at each point. 

“•Jr” a Aj 

J jJ u j, 

3. By means of the values 4 > h'j ' l c , calculating in 
the direction toward the satellite the scattering indicatrix of 
the object and its convolution with the indicatrix of the atmos- 
pheric scattering. 


4. For "pure soil" and the "pure vegetation", calculating by 
means of the direct and the scattered incident radiation, the 
indicatrices and the albedo, the intensity of the radiation escaping 
to the satellite. 


5. Calculating the projective covering at each point from $ ls 
the density of planting, the parameters of the object and the slope 
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of the locality. 


6. After mixing the radiation of the object and the soil 
in proportion to the projective covering, we obtain the radiation 
which escapes from the point to the satellite. 


Then mixing is carried out at points between which the boun- 
dary of objects passes. For example, for the boundary between 



i av,) =■ £(vA' 0; i'U.o-- £ *1 1 

The mixing ratio B is a random variable which is uniform on [0, 


(2/14) 

1 ]. 


Here we take account of the basic phenomena described, for 
example, in [14]. 


Scanning and conversion of radiation into a signal . We shall 
discuss briefly the modeling step under consideration, since this 
question has been thoroughly studied. 


Here the scanning parameters: the diameter of the intake pupil, 
the response and the characteristic signal-to-noise ratio of the 
receiver, the parameters of the preamplifier, etc., serve as input 
data. By the method of finite differences, we solve the equation 
£ for the input signal U(t), where C and R are the 

lit 

load parameters, and i(t) is the input excitation: i(t) = 

+ 1 noise (t ) * * £ [Ti t) 


a 


1 slgnal (t * 

is the receiver response and f u* is its characteristic signal-to 
noise ratio; B is a random Gaussian process with M(B) = 0 corre- 


sponding to the band frequencies _ 


r 


The set of discrete read- 


ings U(t) forms the electrical representation of the frame In 
the given channel. It also serves as the output of the entire 
modeling method. 


Simulation of a probable landscape . This is a rather large 
program unit intended to generate an arrangement of objects, which 


/21 


21 


would not be just a meaningless collection of points, but rather 
would suggest something to the eye. Here objects are chosen which 
have a characteristic form: rivers, reservoirs, settlements, 
forests, fields, etc. In the input data only the desired general 
form of the landscape is specified, and the simulated arrangement 
is random in character, which imparts naturalness to it. 

The second unit is a key factor in satisfying the quality 
criterion for the modeling. Therefore it is realized more care- 
fully and is described in greater detail. The remaining units must 
be regarded as a first approximation. Specialists in the respec- 
tive fields will no doubt improve significantly the quality of 
these units. 

3. THE PROGRAM AND SOME COMPUTATIONAL RESULTS 


The computer program which realizes the proposed method was 
written for the electronic computer YeS-1040. It has a capacity 
of nearly 1000 instructions. The basic computations are performed 
with a frame of dimensions 30 x 30 cells with 10 objects (5 agri- 
cultural crops, a forest, a meadow, roads, water, settlements) and 
with a maximum of 7 parameters (for agricultural crops: the height 
of topography, the soil, the vegetative phase, humidity, produc-' 
tivity, contamination, falling [i.e., of a plant under its own 
weight]). The real relationships between the parameters and the 
albedo, the indicatrix, etc., were taken from [14, 15], and also 
from a collection of journals. Very many of the relationships 
did not turn out successfully. Then we used their heuristic evalu- 
ation. One of the versions of a data bank served as input data. 
Therefore interest was manifested in the volume of data under vari- 
ous conditions. The data bank is presented in Figure 3* The 
largest files involve the relationships of the albedo and the indi- 
catrix in terms of the parameters. 

The five units described in the program are completely inde- 
pendent, and operate consecutively over the files which represent 


/22 


22 



the frame. The computer time of the program is shown in Figure k . 

It was found by extrapolation from the known operating speed of 
each unit. Clearly, to check out the processing procedure it is 
not necessary to use very large frames. In general, the given method 
is designed for laboratory experiments, and the problem of rapidity 
of action is not very critical for it. 

As has been shown, satisfaction of the quality criterion for 
modeling consists of two parts: selection of an e which is suffi- 
ciently small relative to the real data, and sufficiently precise 
reproduction of the input data (the criterion with respect to 3). 
Intrinsically, modeling is characterized by the criterion 3. We 
also checked it out. The input data were chosen with normal distri- 
butions (since they can be characterized by two parameters all told) 
and arbitrary correlation coefficients. The agreement of the simu- 
lated distributions and correlations with the input distributions 
and correlations were vierified respectively by means of the ^ 
and Fisher criteria. It was assumed that they make up part of 
the metric d. The nonstationariness of the' mathematical expectation 
of the input distributions contributed 10% . The results of a 
test with the corresponding significance levels for the deviations 
are presented in the table. 

Clearly, the agreement is good on the whole. This means 
that intrinsically the modeling satisfied the requirements imposed 
on it. In (1.2) a must be chosen on the basis of practical consid- 
erations and an e-network must be constructed which ensures the 
satisfaction of (1.5). 


CONCLUSIONS 


The necessity for mathematically modeling the terrestrial 
surface and the process of remote sensing is shown. The basic prin- 
ciples of the methods of modeling are introduced, in particular, the 
description of objects in terms of their essential parameters. A 
criterion for the quality of modeling is introduced and a practical 
method for satisfying it is proposed. A flexible scheme is suggested /23 
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for modeling a random field with given characteristics, and a 
generalized model of the atmosphere is also given. 

On this basis a full cycle of transformations of solar radiation 
into RS is modeled. A computer program is formulated and tested. 

The testing shows that the modeling procedure complies with the re- 
quirements imposed on it. Moveover, .the following problems are 
posed : 


1. To find a more stringent evalatuion for the quality criterion. 

2. To test on concrete special models the general models for the 
transformations of radiation proposed here. 

3. Practically, using a large number of RS frames, to select 
input data which ensure a good approximation to real frames te- 
net work) . 

4. With the aid of field investigations, to measure all the 
necessary relationships: the dependence of albedo and the indica- 
trix on the parameters, etc. 

5. Finally, to test the processing procedures on synthesized 
and real data and to be convinced with one's own eyes of the effec- 
tiveness of modeling. 

In the present paper these questions are only touched upon due 
to the limited competence of the author or the necessity for ex- 
tensive field investigations. It is to be hoped that if modeling 
interests experimenters, then these problems will be solved with 
expedition. 
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TABLE 


Results of testing the reproduction of input data when modeling the field of parameters (Fixed 
values equal to M +YD' P rovid ed 105K of the readings. The number of readings oscillated between 
50 to 200 for various parameters). 


Designation Input Data 
of parameter 


Modeling Results 


Stationary Cases 


Nonstationary Cases 



! w 


Falling 


A?? : 

0,3 0,003 0,5 t o, 30 
0,008 -0,5' 6,27 
^0~| 0,53 

m n c; 1 


Productivity Q>3 j q^q 

0_,5~ 0,01 ><),5, 0,49 

Contaminationjo ,I| 0,004 0^Z 0,05 

Humidity j o. 2 C>T003 | 0~ 0,18 
Vegetation j 0 ,7|0,0I 0/T 

0,6 0,48 
‘ 0.21 0,27 


) ,02 

0,02 

0762" 

0,015 

0 ,018 

0,009 

0,09 

0,013 

0JDI5 

0,009 

0,007 


0,005 [0,003 

0,0 

pc° jo C03 

oTo io , c 

0, oil 0,003 
0,012 jo, 004 
o,6i i 

0,007 0,003 
0,012 0 . 004 
0,00 3 o', JoT 
6Tol2lo7oo2‘ 


1 6,05 
MM8 
j '-6, 58 
fp.I 
1 

pop 

TfiA 1 

-0,04 

A 55 ! 

jpr 

6,24 


<*(''} 

0,07 

0,05 

0 l 05 

0,09 

0,06 

0,07 

0^03 

0,02 

0,06 

0,04 

0,05" 


Fish- 

Cri- er’s 
terion Cri- 
terion 



15% 7% 

3% 24% 

'1/6 9% 

44# 38# 

11% 26% 


Corre- 

lation 


Error 


Fish- 
er's 
Cri- 
, terion 



- 0,18 0,12 
C,60 0,10 
-0,44 0,11 

-0 ,14 jo ,10 
0,47 0,09 

-0,3$ io.IQ 
0, tO j 0,03 


0 — 

1 

7 — 
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^Dependence of albedo of objects and soil on parameters 
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; ? : results of the operation of the program uriits. The characteristics of the 

a tmo sphere ^urr Tire path of the incident radiation are outlined by dotted lines, and the corresponding 
quantities on the path of the escaping radiation, by dot-and-dash lines. 
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Figure 3. Dependence of the volume of input data V on the number of 
parameters of the ground objects , on the. number of objects N Qb , 

the number of spectral channels N ch , and the degree of polynomial 

used to approximate the different nonlinear functions . On the 

solid lines, N d = 3, and on the dotted lines, = 5* 
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Figure 4. The dependence of computing time t on the dimensions of 

the frames N , N , the number of parameters of the ground objects 
x y 

N , the radius of the correlational neighborhood for the parameter 
P 


R, and the number of spectral channels N c 
mating polynomial N, = 3 • 


The degree of approxi- 


